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The evolution of the gluon plasma produced with saturation initial conditions is calcu- 
i lated via Boltzmann transport theory for nuclear collisions at high energy. The saturation 
scale increases with A and -y/s, and thus we find that the perturbative rescattering rate 
decreases relative to the initial longitudinal expansion rate of the plasma. The effec- 
tive longitudinal pressure remains significantly below the lattice QCD pressure until the 
Q\ . plasma cools to near the confinement scale. Therefore, the transverse energy per unit of 
^ ! rapidity and its dependence on beam energy provides a sensitive test of gluon saturation 
Q I models: the fractional transverse energy loss due to final state interactions is smaller and 
exhibits a weaker energy dependence than if ideal (nondissipative) hydrodynamics applied 
throughout the evolution. 

In collisions of heavy ions at high energy a large number of gluons is liberated from the 
^ , nuclear wave functions. The "plasma" is produced from copious minijet gluons at central 
Q I rapidity, ?/ — 0, with transverse momentum pt > Po- For large po, the produced gluon 
^ ' plasma is dilute. As Pq decreases, however, the density of gluons increases rapidly due to 
the increase of G{x, ~ 2ptI \J s decreases. It has been conjectured |jl],|| that below 

^ . some transverse momentum scale po < Psat the phase-space density of produced gluons 
may saturate since gg ^ g recombination could limit further growth of the structure func- 
tions. Phenomenologically, this condition may arise when gluons (per unit rapidity and 
transverse area) become closely packed and fill the available nuclear interaction transverse 
area. The saturation scale Psat can thus be estimated from dA^(psat) /dy = p'^^^R\/ P, where 
/3 ~ 1. For /5 = 1 the solution reported in EKRT g was psat ~ 0.208A o-^^s where 
Psat and ^/s are in units of GeV, and Ci = dEx / dy / psa,tdN / dy ^ 1.34y4~°-°°'' Our 
focus is to investigate whether the final observed dEj./dy can be used to test the predicted 
A and y/s dependence of the initial dE^p/dy. 

Different gluon saturation models based on classical Yang-Mills equations 0,^ suggest 
that the factor [3 may vary parametrically as /3(psat) = 47rQ;(psat)A'c/c(A''^ — 1), where 
c ~ 1 is a nonperturbative factor proportional to the fraction of the initial gluons in the 
nucleus which are liberated. This factor was recently estimated to be c ~ 1.2 — 1.5. 
On the other hand, PHOBOS and PHENIX data require c ~ 1.9 if one assumes 
entropy conserving hadronization, and for a running with psat; that is, ~ 50% additional 
entropy must be produced during nonideal expansion of the saturated gluon plasma. In 
the approach of 0, the ratio of final to initial multiplicity (or entropy) also grows with 
energy as ~ a'"^^^. This, in turn, suggests that less mechanical work is being performed 
by the longitudinal expansion than predicted by ideal fiuid dynamics 0, and so the Et 
in the final state should be closer to it's initial value. 
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If the expansion proceeds in approximate local equilibrium with pressure p = c^e and 
speed of sound c, then the energy density, e(r), must decrease faster than the expansion 
rate Fexp = 1/r and leads to a bulk transverse energy loss 
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//local equilibrium is maintained during the evolution 6 = In contrast, if the system 
expands too rapidly to maintain local equilibrium, then the effective pressure is reduced 
due to dissipation. The extreme asymptotically free plasma case corresponds to free 
streaming with 5 = 0. Et thus provides an important barometric observable that probes 
the (longitudinal) pressure in the plasma. 

The relevant relaxation rate is given by the fractional energy loss per unit length, 
V^ei = dlogE/dz, which receives a contribution both from elastic and inelastic scattering. 
The relaxation rate is approximately given by |^ 



r,,; ^ 97r«^^ flog - + ^)= K,^97ra^^ log - , (2) 



where we lumped energy loss from radiation into Kin, the inelasticity i^'-factor. Kin = 1 
corresponds to purely elastic scattering, while Kin = 2, for example, corresponds to twice 
the lowest-order elastic scattering rate. 

At the initial time tq = l/psat, s/2 = Cq/po = Cfp^^^. Therefore, noting that the 
comoving gluon density at time tq is po = the ratio of the relaxation rate to the 

expansion rate is given by 

Kin-T—^log- . (3) 
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While Trei oc Psat increases as a power of the energy, the Bjorken boundary conditions force 
the system to expand londitudinally initially also at an increasing rate Te^plro) = Psat- 
The essential quantity that fixes the magnitude of the effective pressure relative to that 
predicted by lattice QCD (LQCD) is the ratio of rates in Eq. (|^), which dimensionally is 
simply a function of a(psat)- The asymptotic freedom property of QCD therefore requires 
that this ratio vanishes as a/s ^ oo. In (j^) the rate of how fast it vanishes is controlled 
by -ft"m«^(Psat)//5(Psat)- Therefore, with saturation initial conditions, asymptotic freedom 
reduces the effective pressure acting at early times r ~ tq and causes the initial evolution 
to deviate from ideal hydrodynamics for a time interval that increases with energy ||^. 

For a more quantitative estimate of the transverse energy loss due to longitudinal 
work, we employ the Boltzmann equation in relaxation time approximation . As initial 
condition, we assume that the initially produced partons have a vanishing longitudinal 
momentum spread in the comoving frame, i.e. we assume a strong correlation between 
space-time rapidity and momentum space rapidity; the accuracy of this approximation 
at the finite RHIC energy remains to be checked within more elaborate models. In any 
case, our gluon plasma thus starts at zero longitudinal pressure, and we shall follow 
its evolution up to the point where it reaches the p/e-ratio from LQCD [§]. In Fig. |l| 
we show the effective longitudinal pressure as a function of the energy density for a/s = 
20, 200, 5400^4 GeV saturation initial conditions [§. Initially p/e starts at zero and remains 
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Figure 1. The ratio of the effective lon- 
gitudinal pressure to the energy density 
along the dynamical path is shown for 
SPS, RHIC, and LHC saturation initial 
conditions^. Solid (dashed) curves are 
for /? = 1 and Ki^, = 1(2). The LQCD 
equation of state [|l^ is also shown for 
comparison. 
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Figure 2. The ratios of the final to the 
initial transverse energy per unit rapid- 
ity are shown as a function of beam en- 
ergy for central Au-|-Au collisions. The 
initial value corresponds to the EKRT 
parametrization {(3 = 1) and the trans- 
port results are for Km = 1,2. The final 
transverse energy for an ideal gas of gluons 
in equilibrium (throughout the evolution) 
is also shown. 



small for a large time relative to tq = 1/psat because the plasma is torn apart by the initial 
rapid longitudinal expansion. The effective pressure approaches the LQCD curve from 
below and reaches it at a time r^^ ^ 1 — 2 fm at RHIC, \/s = 200 A GeV, by which time 
the energy density has dropped by an order of magnitude, ei = ^{tl) = 6.5 — 12 GeV/fm^. 
For LHC y/s = 5400^4 GeV, r^, = 3 — 7 fm during which the energy density falls by almost 
two orders of magnitude to = 9.5 — 21.5 GeV/fm'^. The quoted intervals correspond to 
Kin = 1 — 2 using the EKRT parametrization 0. The contrast between the dynamical 
path followed by the saturated plasma compared to the equilibrium equation of state is 
striking. A qualitatively similar behavior of the early logitudinal pressure has also been 
found from solutions of diffusion equations |TT] . 



The main experimentally observable consequence of the reduced effective pressure is 
shown in the right panel. The ratio Ej,/E^ = Tftf/TQtQ has been obtained from the 
solution of the transport equation assuming free streaming at e/ = 2 GeV/fm^ which 
corresponds roughly to T ~ Tc. tj is determined assuming hydrodynamic expansion from 
the point were the LQCD pressure curve is reached. On the other hand, for an ideal gas 
of gluons produced at time Tq, the final observed transverse energy for 1-1-1 dimensional 
adiabatic expansion would be 



_ -^f^f _ viTfSf-pf) _ Tf 
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Clearly, for Tf ~ one would observe a much smaller transverse energy in the final state 
than in the initial state. Moreover, Ej,/E^ would also have significantly stronger energy 
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dependence such that Ej, deviates more and more from with increasing y/s. In this 
sense isentropic hydrodynamics erases information on the interesting initial conditions via 
this observable. The solutions of the transport equations clearly show a smaller decrease 
of Ej: and of the logarithmic slope, k = dlogii^j^/dlog a/s, due to final state interactions. 
We find that k = 0.50 for the evolution with Kin = 1, k = 0.46 with Kin = 2, while 
K, = 0.40 with isentropic expansion, eq. (^. For comparison, the initial EKRT saturated 
E^^ = TrR^TQeo scales with the higher power k = 0.59. The fractional transverse energy 
loss is thus less dependent on energy than for entropy conserving expansion for which 

f ' 2 

Ej^/E^ oc 1/Tq oc ' . This is due to the increasingly long time spent far from 

equilibrium in Fig. |l| as the beam energy increases. 

The results in Fig. ^ are encouraging from the point of view of searching for evidence 
of gluon saturation in nuclei at high energies. Experimental data on dET/dy or dEx/dr] 
for central Au+Au collisions at RHIC will soon provide a new test of saturation and 
non-saturation models at those energies. Since we predict that dissipative effects reduce 
considerably the effective longitudinal pressure in Fig. |I|, the beam-energy dependence of 
the transverse energy is expected to reflect much more accurately the predicted power law 
dependence of the initial conditions as seen in Fig. 0. We have also shown the preliminary 



result from PHENIX |T|], dEt/dr] ^ 570 GeV for the 2% most central Au+Au events at 
y/s = 130A GeV, scaled to dEt/dy and divided by the initial Et from EKRT. It would be 
most interesting to have data both at higher and lower energy to determine the slope k, 
and compare to the expectation for a saturated plasma. The energy and A systematics 
of the bulk calorimetric observable, dEx/dy, will be a sensitive test of saturation models 
of gluon plasmas produced in the RHIC to LHC energy range. 
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